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Kalman Filtering With Relays Over Wireless Fading Channels 


Alex S. Leong and Daniel E. Quevedo 


Abstract —This note studies the use of relays to improve the perfor¬ 
mance of Kalman filtering over packet dropping links. Packet reception 
prohahilities are governed hy time-varying fading channel gains, and the 
sensor and relay transmit powers. We consider situations with multiple 
sensors and relays, where each relay can either forward one of the 
sensors’ measurements to the gateway/fusion center, or perform a simple 
linear network coding operation on some of the sensor measurements. 
Using an expected error covariance performance measure, we consider 
optimal and suhoptimal methods for finding the hest relay configuration, 
and power control problems for optimizing the Kalman filter perfor¬ 
mance. Our methods show that significant performance gains can be 
obtained through the use of relays, network coding and power control, 
with at least 30-40% less power consumption for a given expected error 
covariance specification. 

I. Introduction 

Since the seminal work of m, which showed the existence of a 
critical threshold on the packet arrival rate for stability of the Kalman 
filter with i.i.d. Bernoulli packet losses, the problem of Kalman 
filtering over unreliable channels has received considerable attention. 
Extensions of include the case of multiple sensors mm , further 
characterizations of the critical threshold (5), probabilistic notions of 
performance (6), the consideration of delays (7), correlated packet 
losses O, O, and power control for stability (9). 

In digital communications, channel coding is often used to improve 
the quality of transmissions over unreliable channels. The concept of 
network coding (To), CD, where in a network with many nodes, 
throughput can be increased by allowing intermediate nodes to 
perform simple operations (such as linear transformations d) on 
its received information, has attracted much attention in recent years. 
Kalman filtering with power control and coding was considered in 
03, CD. The recent work C3l included a study of network coding, 
where one could choose to utilize a relay to perform a network 
coding operation, and the energy tradeoffs involved. The use of relays 
in combating the effects of fading and increasing channel capacity 
has been extensively studied in wireless communications, see e.g. 
d, CD- Indeed, cooperative communications via the use of relays 
has been identified as one of the key enabling technologies for fifth 
generation (5G) mobile networks an. The use of a relay in control 
has been studied in 113 , which showed that for the case of a single 
sensor the stability region for stabilizing an unstable LTI plant can be 
increased in some situations, and also in applications towards control 
of unmanned aerial vehicles (UAVS) Q^. 

In this note, we study remote estimation using relay nodes, and 
investigate what information the individual relays should send to 
the gateway/fusion center. In a related setup considered in 
the relay could only perform network coding that linearly combines 
two of the sensor transmissions using an XOR operation on)- Here, 
we allow for the possibility of the relay combining multiple sensor 
transmissions using XOR operations mi, as well as the possibility 
of the relay forwarding the sensors’ transmissions, which can give 
better performance in certain situations. The current work extends 
our recent contribution (D, which considered the situation with two 
sensors and one relay, and i.i.d. packet dropping links. In this note 
we generalize Q) to multiple sensors and relays, and additionally 
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consider time-varying packet reception probabilities governed by 
time-varying fading channel gains and the transmission power over 
this channel. 

II. System model 

Throughout this note, k will denote the discrete time index, i will 
be used to denote the sensor indices, and I the relay index. The 
process is a discrete time linear system 

Xk+i = Axk+Wk ( 1 ) 

where Xk G K** and is i.i.d. Gaussian with zero mean and 

covariance matrix Q > 0[j The process is observed by M sensors 
with measurements 

yi,k ~ CiXk “t“ Xi^kj i = 1 , . . . , A/ ( 2 ) 

with Di^k G R, Vi and {ui.fe} are i.i.d. Gaussian with zero mean and 
variance Ri > 0, i = The processes {vi^k} and {loji} 

are assumed to be mutually independent, with {A, C) detectable and 
{A, stabilizable, where C = col(Ci,.. ., Cm)- 

We assume that the measurements yi^k have undergone source 
coding and can be grouped into packets of b bits, with each packet 
short enough to be transmitted within one time step. In particular, 
the uniform quantizer of 1201 will be used here. Under the additive 
noise model for quantization (which in general is quite accurate for 
bit rates as low as three bits per sample I21I I. the quantized value of 
yi^k can be written as 

Vi.k ~ yi,k “f Qi,k 

where the quantization noise qi^k has variance with 

j _ 461n2 
“ 3 X 22*’ 

when using the uniform quantizer of 1201 . The measurements are 
transmitted over orthogonal/parallel channels to a gateway, which will 
perform the Kalman filtering operation. Let 'yi^k, for i = 1,..., M, 
be random variables such that yi^k = 1 if is successfully 
transmitted to the gateway by sensor i, and 'yi^k = 0 otherwise. 

Furthermore, there exist L intermediate relay nodes that can be 
used to aid the transmission of the sensor measurements to the 
gateway. Such situations can for instance occur in mesh networks, 
where nodes close to the process will make measurements of the 
process, while the other nodes don’t make measurements but can 
be used to relay the sensor measurements to the gateway M- A 
diagram of the system model for the case of M — 2 sensors and 
L — 2 relays is given in Fig. [D Each relay can listen to a subset of 
the sensor transmissions. Denote I — {1,..., M} as the set of all 
sensors, and X; C T as the set of sensors which relay I can listen to. 
In general the sets Ii,l — 1,..., L will not necessarily be disjoint, 
with possibly multiple relays listening to a given sensor. For i Gli, 
let j. be a random variable such that Cl k ~ ^ if th® transmission at 
time k of sensor i is received by relay I, and Cl k ~ ^ otherwise. The 
relays can perform some simple local processing before transmitting 
over orthogonal channels to the gateway. Let yi^k for 1 = 1,..., L be 
random variables such that yi^k = 1 if transmission at time k from 
relay I to the gateway is successful, and 7 ;,^ = 0 otherwise. 

*We say that a matrix X > 0, if it is positive definite, and X > 0, if it is 
positive semi-definite. 
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Fig. 1. System model for the case of two sensors and two relays 


TABLE I 

Notation for different types oe links 



Channel Gain 

Packet Reception 
Random Variable 

Packet Reception 
Probability 

Sensor i to Gateway 

9i.k 

'li.k 

^i.k 

Relay 1 to Gateway 

9i.k 

ii,k 

\.k 

Sensor i to Relay 1 

h'- 

'^i.k 

^i.k 

Pi.k 


TABLE II 

Boolean expressions for and 02 a . foR two sensors and one 

RELAY, UNDER THREE TYPES OE OPERATIONS 


4>k 

Forward ^ 

Forward 1/2 k 

Send ® yl^ 

01,fc 

7l,fe V(7i,j. AC( j.) 

7l.fc 

7l.feV(7i,s, A72,feACi\,,AC2,fc) 

02,fe 

72,fc 

72,feV(7i,fcAC2,fe) 

72,fe V(7i,s; A7i,fe ACi\fc AC^ j.) 


In this note we will consider a few simple operations that the relay 
can perform|3 A relay can either: a) listen to one of the sensors’ trans¬ 
missions, say sensor i, and forward y1 j, if it is successfully received 
by the relay, or b) listen to a number of the sensors’ transmissions, 
say sensors ii, 12 , • • •, U, and send © ■ ■ ■ © j. if 

Vii k’Vi 2 k’ ■ ■ ■ Vii k been successfully received by the relay, 

where ® is the XOR operation. The XOR operation is commonly 
used in network coding Qni, GD- For instance, if the gateway 
receives both y'^ ,, and y^^, © or both y®and y'^ ,, © i/Jthen 
the gateway can recover both y^ j, and with a binary subtraction. 
In general, given the transmissions received at the gateway, the 
measurements which can be recovered can be determined using 
Gaussian elimination over Z 2 I 3 Determining which sensor/s each 
relay listens to, and which operation each relay uses, is one of the 
key questions addressed in this work, see Section lYl We define a 
relay configuration 

Cf’k = i4’l,k, ■ ■ ■ , 4’L,k) 

at time k as the set of operations (pi^k that each relay uses at time k. 
The set of all possible relay configurations will be denoted by 

The communication channels will be modelled as time-varying 
fading channels |23]. We let gi,k,i = 1,..., M be the channel gains 
at time k from the sensor i to the gateway, gi,k,l ~ 1,... ,L the chan¬ 
nels gains from relay I to the gateway, and £ Xi, I — 1,... ,M 

the channel gains from sensor i to relay 1. We use the block fading 
model HD and assume that {pi.fe}, {gi.k}, {h\ k} vary over time k 
in an i.i.d. manner, with the processes being mutually independent. 
Denote the transmit powers at time k of the sensors and relays by 
Ui.kfi = 1,..., M and ui^k,l ~ 1,...,L respectively. Following 
the model of the packet reception probabilities will depend on 
both the channel gains and transmit powers as follows: We have 
Ai.fe = P(7i.fe = Mgi,k,Ui^k) = p{gi,kUi,k) as the time-varyin£| 
packet reception probabilities from sensor i to the gateway, Ai,fc = 
= llSi.feiWi.fe) = p{gi,kUi^k) the probabilities from relay I 
to the gateway, and = ¥{Clk = = P{h\,kUi,k) the 

probabilities from sensor i to relay 1. Here p{.) : [0, 00 ) —>■ [0,1] is 
a continuous monotonically increasing function whose form depends 
on the particular digital modulation and coding scheme being used 
1231 . For example, in the case of uncoded binary phase shift keying 

2We assume limited computational power at the relays, thus only simple 
operations at a bit-level are considered in this note. If, however, additional 
computational capability is available, then other possibilities include the use 
of more involved network c odin g schemes im or the computation of local 
state estimates at the relays Ea 

^While our use of the XOR operation is similai' to network coding, our 
objectives are not exactly the same. In network coding transmissions are often 
regarded as “successful” only if all packets arrive (eventually) at their intended 
destinations, whereas in our problem even if some packets are not received 
one still can perform state estimation using the available measurements. 

"^This differs from the i.i.d. Bernoulli packet drop model with constant 
packet reception probabilities considered in e.g. (3 


(BPSK) transmission with b bits per packet, p{.) would take the form 



where we assume a packet is successfully received if all b bits 
are succesfully received. However, if there is channel coding and/or 
different digital modulation schemes, p{.) will take on different forms 
(25). In Table U we summarize the notation for the channel gains, 
packet reception random variables and packet reception probabilities 
for the different types of links. 

In addition, there are feedback links from the gateway to the 
sensors and relays which can be used to communicate the relay 
configuration cj)k and power levels m^k and ui^k to be used, see 
Section lYl and Section rviil In this note we will assume that 
transmissions can occur over a much faster time scale than the 
process 0. Thus, delays experienced by the measurements in passing 
through intermediate relay nodes will be ignored. For instance, in 
the industrial wireless sensor networks standard WirelessHART (261 . 
transmissions between nodes would typically take around 10 ms, 
whereas for many estimation and control applications the process 
time constant might be 1 sec or more. 

III. Kalman filter with packet drops and relays 

Let 9i^k,i = 1, ... ,M be random variables such that Oi^k = 1 
if t/?j, can be reconstructed at the gateway, and Oi^k = 0 otherwise. 
Note that in general Oi^k and Oj^k, i j are not independent. Values 
of Oi^k for different relay configurations and combinations of yi^k, 
li,k, Ci k^ can be written in Boolean algebra form. For example, in 
Table [ft] we give the Boolean expressions for 9i^k and 02,fc in the 
case of two sensors and one relay, where we use the notation A to 
denote logical ‘and’ and V to denote logical ‘or’. Now define 

C’fc COl(0l^fcC’l,..., 9 ^ kC M ), yk • • • 5 ^ .ky ^ k ) 5 

■ ■ ■ 7 yk., Coj . . . , Ck\ 

Pk+l\k (Xfe-|_l ^ky-l\k^ 11/0; • • • i Pk j C’Oi • • ■ 5 C’fc] ■ 

(4) 

The associated Kalman filter equations which are run at the gateway 
can be written as (see e.g. HU) 

— 1 © XCkiyk CklCk\k—l') 

Pk + l\k = ^Pklk-lA^ + Q — KkCkPk\k-l^ 

where Kk = APk\k-iCl{CkPk\k-iCl +Rk)~^, with Rk = 
diag(.Ri,fe,..., Ru.k) = diag(7?i+5(,E[i/i j.],..., RM+5b^[y'ii,k\)^ 
similar to (H In the sequel, we will also call Pk = Pk\k-i- 

Remark III.l. As in idl, the Kalman filter 0 uses all successfully 
reconstructed measurements, but those measurements where 9i^k = 0 
are not taken into account. 
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IV. Performance oe the Kalman filter with relays 
In Section |n] we have proposed that each relay can either listen 
to transmissions from one of the sensors which it then forwards 
to the gateway, or it listens to a number of sensors and performs 
an XOR operation that is then sent to the gateway. We wish to 
investigate which operation each relay should use, and which sensors 
each relay should listen to, i.e. determining the relay configuration 
in order to give the best performance for the Kalman filter. 
This section presents some preliminary results on the performance 
of the Kalman filter, with optimal relay configuration selection to be 
studied in Section lYl For time-varying Kalman filters with packet 
drops, one commonly used performance measure is the expected 
error covariance. In this note we consider the problem of optimal 
relay configuration selection in order to minimize the one step ahead 
expected error covariance E[Pfc+i|Pfe, gj,, </)fc], where 

gfc — {si.fei • ■ • I gM,k,gi,k, ■ ■ ■ ,gL,k, fi-i.fc, • • •, hM,k} (6) 


represents the channel gains at time k, which in turn will determine 
the packet reception probabilities Ai,fe, \i,k, Pi^k^i = li ■ • ■ i M, I = 
1,..., L. In order to compute E[Pfc+i|Pfc, gf., (f)k], we will further 
assume that full channel state information (CSI) at the receiver is 
available, so that g^ is known at the gateway|3 Define 


fk{X)^AXA'^+Q-E 


AXCl {CkXCl+Rk) ^CkXA^\g^,<i}k 

(7)' 


where the expectation is with respect to 6i^k, ■ ■ ■ ,dM,k in the 
definition of Ck in @. Equivalently, we can write fk{X) as: 

fk{X) = AXA'^+Q- ^E [n n (1 - I Sic ’ 

ICI i€l j^I 

X (^ C { I ) XC { lf + Rk { I ) yC { I ) XA '^ 

( 8 ) 


Cl,k 

_A 

{Ai.fc, . 

. . , XM,k., 

, Ai,fe, . 

• • ? ^l,k ; • • • ? ^L,k 

1 

? Pl,k, • 

L 

■ ■ , pM,kS 

Idl^k 

_A 


. . , XM,k. 

7 M,k, . 

• • 5 l^l,k ? ■ ■ • ? ^L,k 

1 

5 Pl,kf ■ 

L 

■ ■ 7 pM,k S 

Ti\,k 

_A 


. ■ , XM,k 

) , ■ 

■ ■ , ^L,k^Pl,k^ ■ ■ ■ 

/ 

? Pi,k ? • 

L Y 
• • 7 PM,k S 

Si,k 

_A 

{Ai,fe, . 

. ■ , XM,k 

I ■ 

■ ■ , ^L,kj Pl^kt ■ ■ ■ 

1 

^ ^i^ki • 

L -\ 

■ ■ 7 PM,kS 

Lemma 

IV.l. Let fx,^,{ 

.) be defined by fk{.) in 

((T] when the links 


have packet reception probabilities Xi^k- Then, irrespective of which 
relay configuration is used, Vi = 1,..., M, V/ = 1,..., L, 

Ai,fe < pi,k => ^ 

\i,k < fii,k ^ 

Pi,k < (^i,k ^ f-Tp, ^ ^ fs[ j, (^) 

Proof: Consider the case Xi^k < Pi,k- Recall that Bernoulli 
random variables can be generated from P(0,1) uniform random 
variables, by comparing the uniform random variable with the 
probability that the Bernoulli random variable is equal to one, i.e.. 
7i,fc = 1 when u < Xi,k, and = 0 otherwise, where u is 
P(0,1). Let LU denote an outcome corresponding to N indepen¬ 
dent realizations of P(0,1) random variables, where N is equal 
to the total number of packet dropping links. For each lo, one 
can generate corresponding independent Bernoulli random variables 
'n,k,---,yM,k,yi,k,---,yL,k,(lk,---,CM,k- One can then con¬ 
struct the Bernoulli random variables 8i,k, ■ ■ ■ ,8M,k, and hence Ck 
as in 1 ( 4 ). 

Let C'c- {uj) be the matrix Ck when using packet reception 
probabilities Ci^k and Cui be the matrix Ck when using packet 
reception probabilities Ui,k. Now note that if Oj^k{io) = 1 using the 
packet reception probabilities jCi,k, then we also have 9j^k{oj) ~ 1 
when using the packet reception probabilities Ui,k, by the way in 
which 9j^k{oj) is constructed, and since an increase in the packet 
reception probability of any link cannot decrease the probability of 
reconstructing any of the sensor measurements. Hence 


where C{I) = col({Ci,i € /}), Rkil) = diag{{Ri^k,i £ 7}), and 
Oi.fc, i = 1, ■ ■ ■, M are random variables with the same distributions 
as The quantities E ©i.fc can 

be computed in terms of the packet reception probabilities Xi,k, Xi,k, 
Pi^kA = li ■ ■ ■: XI, I = 1,..., L. A systematic procedure for doing 
this is as follows: 

1) Write out the Boolean expression 

/\8i,k /\{-^0j^k) (9) 

iei j^i 

where each Oi^k is written as a Boolean expression, -^9j^k denotes 
the negation of the Boolean expression for Oj^k, and the notation 
Aie/ Oi,k = Oij^,k A ■ ■ ■ A 6i„^k A ... for indices i„ £ I. 

2) Convert the Boolean expression © into the sum of products 
normal form (27]. (Note that this can be done systematically.) 

3) E[nie/“ ^j,k)\gk,(f>k]is then given by taking 

the sum of products normal form of (3, and replacing A with 
multiplication, V with addition, with Xi^k, ~'yi,k with 1 — Xi^k, 
yi^k with Xi,k, -•fi.k with 1 - Xi,k, CU and ^ with 

1 — Pi,k‘ 

By step 2), each term in the sum will correspond to a dis¬ 
tinct entry of the truth table for ..., thus allowing 

® [riig/ ©i.fe - 0fA)|gfci to be easily calculated. 

We now give a result on how the packet reception probabilities 
affect the expected error covariance E[Pk+i\Pk,gk^*t^k\ = fk{Pk)- 
First denote: 


AXCc,,k{i^f {Cc,,k{i^)XCc,,k{^f+Rk) Cc,,kHXA'^ 

> AXCu.Jujf yu^ ,{u})XCu,Ju}f+Rk)~Cu,Juj)XA'^ 

( 10 ) 

Since l llOb holds for all uj, we have 


E 


AX{Cc,,kf {Cc,,kX{Cc,,kf+Rk) Cc^.XAy^fik 


> E 


AX{Cu,,,f {Cu,,kX{Cu,,kf+RkyCu,,kXA'^\gk,cf,k 


which shows that fc^ ^ (X) > fu^ ^ (-V). The other two cases can be 
proved in a similar manner. ■ 

Lemma IIV. 1 1 says that increasing the packet arrival rate on any 
one of the packet dropping links will give a decrease in the value of 
the one step ahead expected error covariance, no matter which relay 
configuration is used. 

We also wish to determine how the system parameters such as 
C\,, Cm,Ri,..., Rm will affect the values of the expected error 
covariance. Denote: 


Ai = {Cl,..., Ci,..., Cm, Ri, ■ ■., Ri,..., Rm} 

Bi = (Cl ,..., Di,..., Cm, Ri, ..., Si,..., Rm] 

Lemma IV.2. Let fy.^k{X) be defined by fk{.) in ([S] when the 
sensors have parameters yi. Then, irrespective of which relay con¬ 
figuration is used, Vi = 1,..., M, 


T'i,k — {Al,fc,...,Ai,A:,..., XM,k, Al,fe, • • . , Xc.k, Pl.fc 5 • • • ; PM,k} 
Idi^k {Al.fc, ■ ■ ■ , pi,k , ■ ■ ■ , XM,kj Ai^/c, . . . , Xc.k^ Pl,fc 7 • • • 5 pM,k} 

^In practice this can be achiev ed u sing channel estimation and prediction 
algorithms, see references in di, na 


CjRfiXCi < DfS-^kDi => fMX) > fB,{X) 

The proof of Lemma liy.2l is omitted; a related result is proved 
in (T|. The quantity Cf R~lCi can be regarded as the signal-to- 
noise ratio (SNR) for sensor i (which includes the quantization 
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noise contribution). Thus Lemma IIV.21 shows that the expected error 
covariance is monotonic in the SNRs Cf R~^Ci. 

V. Optimal relay configuration selection 

We now wish to address the question of determining which config¬ 
urations for the relays will give the best Kalman filter performance. 
Suppose for now that the sensor transmit powers Ui^k,i = 1, ■ ■ ■, M 
and relay transmit powers are given or fixed (The 

more difficult problem of jointly optimizing the relay configuration 
and transmission powers will be considered in Section rvlll ). We 
wish to choose at each time instant k, the relay configuration (pl. 
that minimizes TrE[Pfc+i|Pfc, gj., cjik], i.e. 

</>fc = argmin TrE[Pfc+i|Pfc, g^., </>*:]. (11) 

4’k(Pk,t.k)<^^ 

where E[Pfe+i|Pfc, g^., </)fc] = fk{Pk), see ([U. 

A. Optimal relay configuration selection 

Problem GB can, in principle, be solved by exhaustive search at 
the gateway. The optimal configuration can then be fed back to the 
relays. We will characterize the number of relay configurations that 
need to be checked at each time instant for exhaustive search. 

Lemma V.l. Let Ii be the set of sensors that relay I can listen to, 
and let Mi = \Ti \ denote the cardinality ofli. Suppose that there are 
no restrictions on how many relays listen to the same sensor. Then 
the number of possible relay configurations for <fik is 

1^1=n (2"^' -1) (12) 

1=1 

Proof: First fix a relay I, which can listen to Mi of the sensors. 
This relay can either forward any one of the sensor transmissions, 
or perform the XOR operation on two or more of the sensor 
transmissions it listens to, resulting in Mi + {^‘) + (^‘) + • • • + 
(m*) ~ 2**^* — 1 possible operations. If there are no restrictions on 
multiple relays listening to the same sensor, then by the multiplication 
principle the number of relay configurations is HrLi “ l)- * 
We thus see that the number of configurations that needs to be 
checked is in the worst case (where each relay can listen to all 
sensors) exponential in the number of relays L and number of sensors 
M. However, in practice, due to geographical considerations, the 
number of sensors Mi that each sensor I listens to is often small, 
e.g. in l28l . 1291 it is assumed that Mi < 3. 

B. Stability of Kalman filtering with relay configuration selection 

We now wish to give a condition for stability of the Kalman filter 
with optimal relay configuration selection. 

Definition 1. The Kalman filter is said to be exponentially bounded 
if there exist finite constants a and and p € [0, 1), such that 
TrE[Pfc] < ap'^ + Vfc. 

Theorem V.2. Let {sfc} be a stochastic process such that Sk = ^ if 
Ck is full rank, and Sk =0 otherwise. Suppose there exists a policy 
4>\gk) dependent only on g^, such that 

\\A\\^E[nsk = 0\gk,ct>Hgk))] < h (13) 

where | |4l| | is the spectral norm of A. Then the Kalman filter using the 
configurations obtained from problem {77]) is exponentially bounded. 

Proof: Since the distribution of Ck depends on Pk and 
gj., and gj, is independent in time and of Pk, we have 
E(C'fe|Pfe,Pfe-i,---,Po) = P(C'fe|Pfe). Then^ by {3, the pro¬ 
cess {Pfe} is Markovian. Now define 14 = TrP*;. Since in 


the relay configuration selection problem (E) we are minimizing 
E{Vk+i\Pk,gk,<l>k{Pk,gk)}, we have 

E{Vk+i\Pk}=E[E{Vk+i\Pk,gk,tl>UPk,gk)}] 

<E[E{Vk+i\Pk,gk,cf>\gk)}] 

= E[E{Vk+i\Pk,gk,cl>\gk),Sk = l}P{s. = l\Pk,gk,tl>Hgk)} 

+ E{Vk+i\Pk,gk,(pHgk),Sk = 0}P{sfe = 0\Pk,gk,(l>Hgk)}] 

< IK + (||4||Vfe+TrQ)E[P{s, = 0|g„</.“(g,)}] 

where the last inequality is shown using similar arguments to 
1^ . 1131 . and IK is a positive constant. If | |Al| |^E[P(sfe = 
OiSfci < 1 we may then use a stochastic Lyapunov function 

argument similar to ITtI to show that E{Kfe|Po} < ap^ + P,'ik for 
some p £ [0,1) and constants a and /3, which establishes exponential 
boundedness of the Kalman filter. ■ 

Theorem I V.2I thus provides a sufficient condition for Kalman filter 
stability dependent on the system matrix A and the distributions of 
the channel gains g,.. 

Example V.3. Consider the case of two sensors and one relay, with 
Ck being full rank only when 6i^k = 02,k = 1. Then 

P(^l,fe = 1, 02,k = l|gfc, <pk) = E[0l,fe02,fe|gj., Cfk] 

= M,k^2,k + (1 — M,k)M,kM,kpl,k 

Suppose we choose cj)^ to be the suboptimal policy that always 
forwards t/j j., and with the transmit powers ui,k = U 2 ,k = ui,fc = 1- 
The condition {71 then becomes 

E[P(sfc = 0|gfc, = j (i-p{gi,k)pig2,k) 

-('^-p{gi,k))p{g2,k)p{gi,k)p{hl^k)'^d¥{g^) < 

which can be checked by numerically computing the integral for 
specific functions p{.) and probability distributions P(gfc). If con¬ 
dition {21 is satisfied for this suboptimal policy, then by Theorem 
\V.2\ the Kalman filter using the optimal relay configurations will also 
be stable. 

C. Suboptimal relay configuration selection 

Lemma rvTI has shown that the optimal way of choosing the relay 
configuration by checking each configuration is exponential in the 
number of relays L, which is computationally intensive when L is 
large. To reduce computational complexity, a suboptimal method for 
determining a relay configuration is to optimize the operation of each 
relay I independently of each other. The motivation for this method 
is that sometimes other relays may become unavailable, thus one 
should optimize the performance of each relay irrespective of what 
the other relays are doing. Specifically, consider subsets p CP;. Let 
C{Ii) = col({Ci,i G P}), Pfc(p) = diag({Pi,fc,i G P}), and 

fpkiX)^AXA^+Q- 

qci, i€Ii j^Ii 

X AXCiIif(^CiIi)XC{Iif+RkiIi)Jc{Ii)XA'^ 

where the terms E* (1 “ 0j.fc)|gfei <('<,*:] are com¬ 

puted assuming that relay I is the only relay available. One then 
computes fi,k{Pk) for each of the operations fipk that relay I can 
perform, with the one that gives the smallest value of Ti{fi,k{Pk)) 
then chosen. This optimization can be carried out for each relay 
independently of the other relays. The number of configurations that 
need to be checked at each time step k is then ~ l)> 

which (compare to ) ll2b ) is no longer exponential in the number of 
relays P, and with Mi usually being small as mentioned at the end 
of Section IV-AI 
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VI. A SPECIAL CASE 

Here we consider a special case where additional analytical results 
can be derived. We will study the effects of varying the packet 
reception probabilities and signal-to-noise ratios, which will provide 
some insight into the general behaviour. 

Recall that in the case of two sensors and one relay, the relay 
can either i) forward if it is received, ii) forward y 2 ,k if it is 
received, or iii) send ?/i © 2/2 fc Vi k k received. 

We will consider a scalar example with = p\ k = P 2 k ~ Vfc, 
corresponding to the case where the channels from the relay to the 
gateway, and from the sensors to the relay, are error free. We will call 
A = a, Q ^ q. Cl = Cl, C 2 = C 2 , Ri = ri, R 2 = r 2 , Ri,k = ci.fc, 
R 2 ,k = r 2 ,k- We have 




2 7^ -^2 (1 — -^2 r ^ Q 

a^Pk + q-— -- „ , ’ 1 ^ forward 

4/c,k+4/-2,k 

2 n -^1 (1 — -^1 c 1 0 

a^Pk^q- p , ’ 1 ^ ;- - ^ forward 


Pk Pq~ 


(■^1 fc — k^'2 

p I ,, ^ 


, send yl i^i 


We want to see under what conditions the XOR operation j, © 
1/2 f. outperforms forwarding of measurements. First let the terms 
Ci/ci,fe and c|/r 2 ,fe, which can be regarded as the signal-to-noise 
ratios (SNRs) of sensors 1 and 2, be fixed. From (O we see that 
sending yl k® V^k better than forwarding yf ^ when 


have been observed in simulations for more general cases of packet 
dropping links to and from the relays, and with larger networks. 
Similar behaviour has also been reported in the network coding 
literature f2^ . 

VII. Relay coneiguration selection and power control 

In Section|V]the sensor and relay transmit powers were assumed to 
be fixed. However, the presence of time-varying fading channels will 
also allow for the use of power control techniques to further improve 
performance. In this section we present one possible formulation 
which optimizes the estimation performance subject to a sum of 
transmit powers constraint. 

As in Section m we assume that full channel state information 
(CSI) is available at the receiver, with gj, in representing the 
set of all channel gains at time k. The transmit powers of the 
sensors and relays can then depend on both the instantaneous 
channel gains gj, and the error covariance Pk, with these 
transmit powers being computed at the gateway (which is 
assumed to have more computational resources than the sensors 
and relays) and fed back to the sensors and relays before 
transmission occurs at the next time step. Denote Uk{gk,Pk) — 
{Ul,k{gk^ Pk), ■ ■ ■ , UM,k{gk,Pk), Ul,k{gk,Pk), ■■ ■ , UL,k{gk, Pk)} 
as the set of all transmit powers at time k. 

A. Optimal power control for a given relay configuration 

For a given relay configuration, we pose the following power 
control problem: 


P I _ 1 

Ai,fe(l — A2,fe) 1 — A2,fe \ 

A.. ■ 


P _j_ _ 1 

<^'iPl,k^c\lr2,k 

Similarly, sending Hi k® vtk better than forwarding y'^ 


ii/c,k 

(1 

wher 


A2,fc > 


Pk + 


Cl/t‘l,lc+ci/'r2,fc 


Pk + 


ip2,k 


(16) 


From (dS-ldlll, sending yl k®y 2 k best when both packet reception 
probabilities Ai,fc and A 2 ,fc are above certain thresholds, which in turn 
implies that the instantaneous channel gains gi^k and g 2 ,k need to be 
above some thresholds. Thus, for lower quality channels, forwarding 
of measurements gives better performance than the network coding 
operation. The intuitive explanation for this is that when the gateway 
receives Vi k® V 2 k, b nssds one other measurement (either yf f, or 
j.) in order to be useful. In contrast, if the relay forwards yf or 
j., this value (if received at the gateway) can be used even if the 
direct transmissions from the sensors are lost. 

Alternatively, regard Ai.fc and X 2 ,k as being fixed. Rewriting (BJ 
and OUl, sending Vi k® k better than forwarding y\ when 


Ai,fe 

Ci/ri.fc 


1 

C?/Tl,fc + Cllf 2 ,k 


> Pk{'^ 


Ai,fc), 


and sending yl k® vtk better than forwarding 1/2 k when 


(17) 


X2,k 

4/r2,k 


_ 1 _ 

C?/Tl,fc + ci/T2,fc 


> Pki^ 


X2,k) 


(18) 


For fixed Ai,^; and X 2 ,k, we see that if either Cijfi^k or c|/r 2 ,fe 
is sufficiently large, then conditions 01} and Gl cannot both be 
simultaneously satisfied. The intuitive reason for requiring the signal- 
to-noise ratios to be small in order for network coding to give benefits, 
is that the relative performance gains by having both measurements 
available at the gateway (vs. just one of the measurements) is greater 
at low SNRs than high SNRs. 

These qualitative observations, that the XOR operation needs 
channel conditions to be good, or for the signal-to-noise ratios to 
be low, in order to give benefits over forwarding of measurements. 


min TrE[Pfc+i|Pfc,gj,,(/)fc] 

'ikiikPk) 

M L (19) 

S.t. ^ Ui ^ k {% k , Pk ) +^ Ul , k ( gk , Pk ) < Wtot 
i=l i=l 

which minimizes the expected one-step ahead error covariance subject 
to the sum power Y)T^i Ui,k{gk,Pk) + YH^i ui,k{gk,Pk) being less 
than Utot- Due to the objective being a complicated nonlinear function 
of the transmit powers Ufc, the optimization problem 01 is in general 
non-convex and will need to be solved numerically using global 
optimization algorithms 1301 . 


B. Joint relay configuration selection and power control 

Problem 01 is for a given relay configuration. To optimally 
choose both the relay configuration and transmission powers, we can 
in principle solve HrLi ~ l) instances of problem 01 at each 
time step (for each of the configurations, see Lemma lyTt . and choose 
the relay configuration that gives the smallest value for the objective 
function, which however is very computationally intensive. 

A less computationally intensive suboptimal scheme is to first 
choose a relay configuration by assuming a simple power allocation 
(e.g. that the total power Uat is equally divided between the sensors 
and relays), and using the suboptimal method of Section IV-CI to 
choose a relay configuration. For this chosen relay configuration, we 
then further optimize the transmission powers by solving the power 
control problem 01 - 

VHI. Numerical studies 

We first look at the performance differences between the optimal 
relay configuration selection and the suboptimal methods of Section 
13 We consider a situation with two sensors and two relays, where 
each of the relays can listen to both sensor transmissions, see Fig. 
[T] We consider the scalar case with a = 0.95, g = 1, ci = C 2 = 1, 
ri — r 2 = 1. For simplicity, we assume that the links from the 
sensors to the relays are perfect, with the fading channels (from 
the sensors to gateway, and from the relays to gatewary) being 
exponentially distributed with mean 1, which models the case of 
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average sum power 

Fig. 2. Optimal and suboptimal relay configuration selection 


Rayleigh fading (2^. Similar to (3, we assume that the digital 
communication uses BPSK transmission with 6 = 6 bits per packet, 
so that the function p{.) in Section|n]has the form (O. We distribute 
the transmit powers equally between the sensors and relays. Fig. 
plots the average sum power and expected error covariance E[Pfc] 
(obtained by time averaging {xk - £k\k-i)ixk - Xk\k-i)'^ over 
10000 Monte Carlo iterations), for the optimal and suboptimal relay 
configuration selection methods. For comparison we also plot the 
performance for the cases of: 1) no relay, 2) a scheme where the 
relay always performs the XOR operation (m, 3) a scheme where 
the relay sends the two sensor measurements with less accuracy 
by removing half the bitfl and 4) a scheme where the gateway 
can ask for each lost transmission to be retransmitted oncfl In 
each case, the expected error covariance decreases as the average 
power is increase. Since by (0 larger powers imply higher packet 
reception probabilities, this behaviour is in agreement with Lemma 
IIV. 1 1 We also see that the suboptimal method that optimizes each 
relay separately gives very close performance to the optimal method, 
and significantly outperforms the other schemes. 

We next consider the case of two sensors and one relay, with 
Rayleigh fading for each of the fading channels. We choose gi,fc, g 2 ,k 
to have mean 1, while gi,k, h\^k, ( 12 ,fc have mean 4. This models the 
case where power decays in free space as 1/d^ where d is the distance 
from the transmitter |23l, with the relay located approximately 
halfway between the sensors and gateway. In Fig. [3] we plot E[Pfe] 
(obtained by time averaging {xk - Xk\k-i)ixk - over 

10000 Monte Carlo iterations) for different sum powers, obtained 
by solving problem d using the fmincon routine in Matlab for 
each relay configuration and selecting the best one. We also plot the 
performance of the suboptimal scheme where a relay configuration 
is first chosen (assuming equal power allocation) and then power 
control is performed. We compare this with the case where there is 
no power control, with the sensors and relay using the same transmit 
power at all times, but with the best relay configuration chosen at 
each time step. Additionally, we plot the case where the relay always 
performs the XOR operation, and the case without a relay but with 
power control. We see that doing power control gives significant 
benefits, with the best performance when one optimizes both the 
relay configuration and transmit powers. The suboptimal scheme 

®Here the relay removes from each quantized sensor measurement the least 
significant 6/2 bits. Then in the Kalman filter equations l(5}, for the case 
where the direct transmissions are not successful, in the expression for Rk 
we replace S), = with 5(,/2 = ^^^ 2 ^ whenever we use a measurement 

where half the bits have been removed. 

^Here we assume that additional transmit power (same as the power for 
a single transmission) is used for each retransmission, with a successfully 
retransmitted measurement (from time k) available to the Kalman filter at 
time A: + 1, which now utilizes a buffer similar to 0. 



average sum power 


Fig. 3. Power control and relay configuration selection 


where a relay configuration is first chosen by assuming equal power 
allocation, and then the powers are optimized, performs very close to 
the optimal scheme. Comparing the plots using power control, we see 
that for a given expected error covariance the average power required 
is signficantly less (at least 30-40%) when a relay is used. 

IX. Conclusion 

This note has studied the use of relays for Kalman filtering with 
multiple sensors over packet dropping links, where the packet recep¬ 
tion probabilities are governed by fading channel gains and sensor 
and relay transmit powers. By allowing relays to either forward one of 
the sensor’s measurements or perform a network coding operation, 
we have considered the problem of determining the optimal relay 
configuration at each time step, together with a simpler suboptimal 
method. We have also studied the use of power control in addition to 
selecting the best relay configuration, to further improve performance. 
Numerical results have demonstrated that the use of relays can lead 
to power savings of at least 30-40%. Future work will include studies 
of other strategies that relays can perform, and extending our setup 
to utilize relays in multi-hop networks [m. 
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